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FUEL CONSERVATIVE GUIDANCE CONCEPT FOR SHIPBOARD LANDING 


OF POWERED-LIFT AIRCRAFT 

David N. Warner, Jr., Leonard A. McGee, John D. McLean,* and Gregory K. Schmidt* 

Ames Research Center 


SUMMARY 


A simulation study was undertaken to investigate the application of Energy Con- 
servative Guidance (EGG) software, developed at NASA Ames Research Center, to improve 
the time and fuel efficiency of powered-lift airplanes operating from aircraft 
carriers at sea. The ECG software system consists of a set of algorithms whose coef- 
ficients and parameter limits are those of the Quiet Short-Haul Research Aircraft 
(QSRA) . When a flightpath is indicated by a set of initial conditions for the air- 
craft and a set of positional waypoints with associated airspeeds, the ECG software 
synthesizes the necessary guidance commands to optimize fuel and time along the 
specified path. A major feature of the ECG system is the ability to synthesize a 
trajectory that will allow the aircraft to capture the specified path at any waypoint 
with the desired heading and airspeed from an arbitrary set of initial conditions. In 
this study, five paths were identified and studied. The first path closely follows 
the manual approach procedures specified in the U.S. Navy CV NATOPS MANUAL for prop 
and turboprop aircraft using tactical air navigation aid (TACAN) as the major area 
navigation aid. Each of the four remaining paths were established to successively 
remove the manual restrictions from the path. These paths demonstrate the ECG sys- 
tem’s ability to save flight time and fuel by more efficiently managing the aircraft’s 
capabilities. 

Results of this simulation study show that when restrictions on the approach 
flightpath imposed for manual operation are removed completely, fuel consumption 
during the approach was reduced by as much as 49% (610 lb fuel) and the time required 
to fly the flightpath was reduced by as much as 41% (5 min) . When it is possible to 
remove only a portion of the operational constraints, somewhat lesser, but still sig- 
nificant, savings in both fuel and time were realized using the ECG synthesis soft- 
ware. Savings due to ECG were produced by (1) shortening the total flight time; 

(2) keeping the airspeed high as long as possible to minimize time spent flying in a 
regime in which more engine thrust is required for lift to aid the aerodynamic lift; 

(3) minimizing time spent flying at constant altitude at slow airspeeds; and (4) syn- 
thesizing a path from any location for a direct approach to landing without entering 
a holding pattern or other fixed approach path. 


INTRODUCTION 


Aircraft operating from ships at sea must constantly monitor fuel consumption to 
ensure a safe return to the ship and to accomplish the prescribed mission at the 
least cost. It has been demonstrated that energy conservation techniques managed by 
computer-driven displays to the pilot or by completely automatic systems are capable 
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of substantial fuel and time savings, when operational constraints are minimized and 
full advantage is taken of the aircraft capabilities. In situations in which only a 
limited number of the operational constraints can be removed, the savings are also 
limited, but when summed over a large number of flights, such as in aircraft carrier 
operations, the fuel savings can still be substantial over a period of time. 

Recent studies at NASA Ames Research Center, described in references 1, 2, and 3, 
have shown the effectiveness of an on-board software system designed to achieve 
energy conservation. The Energy Conservative Guidance (ECG) software system, which 
consists of a set of algorithms designed for the NASA Ames Quiet Short-Haul Research 
Aircraft (QSRA) , synthesizes flightpaths composed of capture and reference segments 
that optimize fuel use along a route described by fixed waypoints. A !, command table" 
is generated and stored in the computer which calculates a schedule of settings of the 
outboard flaps, the upper-surface-blown (USB) flaps, and throttle at each waypoint and 
at selected intermediate points. This information could be displayed to the pilot for 
manual actions at appropriate times or could be used by a fully automatic guidance 
system, thereby allowing the approach to be flown by the on-board computer system. 

The time to fly the synthesized flightpath is determined during the synthesis proce- 
dure and may be used for orderly and fuel-efficient control of traffic by the air 
controllers on the carrier. 

This paper presents a simulation study of the QSRA aircraft making energy- 
conservative landings on an aircraft carrier deck. The feasibility of landing the 
QSRA on the carrier USS Kitty Hawk has been demonstrated (ref. 4). The carrier is 
under way at a constant speed and direction while heading its canted deck into the 
wind. The fuel and time used by the aircraft are computed for five different 
approaches, each of which is defined by a set of waypoints. In the first two 
approaches, the aircraft enters a standard traffic pattern, which is translating 
along with the ship. This pattern is described by a set of waypoints in the computer 
aboard the aircraft and by the initial conditions of the aircraft at the time the 
flightpath synthesis begins. From these initial conditions and fixed waypoints, the 
computer's ECG software synthesizes a reference flightpath. This flightpath consists 
of a capture path from the aircraft’s current position and heading to a fixed path 
defined by the given waypoints. Time histories of various parameters are generated as 
are the elapsed time and the distance flown. The second flightpath was the same as 
the first except that the speed profile following the turn which aligns the aircraft 
with the centerline of the canted deck is flown at a higher speed. In the third 
flightpath, only two waypoints which describe the final straight-in portion of the 
approach are specified. The first Is at 3 distance-measuring equipment (DME) n. mi., 
and the second is at touchdown. The fourth flightpath is the same as third except 
that the first waypoint is at 1 DME n. mi. The fifth flightpath differed from the 
fourth in that the aircraft's initial X-position is the negative of the previous 
initial X-positions, the heading was -90°, and the first waypoint was at 
0.5 DME n. mi. from touchdown. Thus, each successive flightpath allowed the ECG 
software in the aircraft computer to have greater freedom to utilize the aircraft’s 
aerodynamic capabilities by specifying flightpaths with successively longer capture 
paths and shorter fixed paths in order to synthesize reference paths offering improved 
fuel conservation as well as shorter flight times. 

Navigation information is assumed to be initially supplied by a tactical air 
navigation air (TACAN) system aboard the aircraft carrier. This information would be 
used throughout the approach until visual contact with the landing deck is made by the 
pilot, whereupon the pilot would manually make the landing on the deck. 
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This study compares the fuel used and time required far a landing on the carrier 
for each of the five approach paths. Comparison of results will demonstrate that ECG 
software in an on-board computer can be effective in saving fuel in an aircraft 
carrier operational environment and that fuel savings increase substantially as oper- 
ational restrictions are removed from the flightpath. 


ENERGY CONSERVATIVE GUIDANCE SYSTEM 


The ECG system is briefly described here. A detailed description is given in 
appendices A, B, and C. The purpose of the ECG software is to synthesize a flyable 
reference path from the aircraft’s current position and course to a desired final 
position and course. The synthesized horizontal path consists of a fixed path speci- 
fied by a set of input waypoints and a capture path from the aircraft’s current 
position and course to the position and course at the capture waypoint on the fixed 
path. The capture path consists nominally of an initial turn, a straight segment, 
and a final turn. In special cases, one or both turns may become a straight line. 

In this study, the waypoints specifying the reference path are defined relative to a 
moving-ship-centered coordinate system. 

In an operational situation, input data would be supplied by the pilot, or the 
on-board navigation system, or both. Furthermore, the pilot would have the option at 
any time to change the input variables and request the ECG system to synthesize a new 
reference path which would meet modified operational criteria. The pertinent input 
data used to specify the path’s, the aircraft’s, and the ship’s initial conditions 
are described as follows; 

X(I) ,Y(I) ,Z(I) cartesian position for the Ith waypoint with respect to 

the ship 

R(I) turn radius at the Ith waypoint 

VNOM(I) ,VMAX(I) ,VMIN (I) nominal, maximum, and minimum indicated airspeed at the Ith 

waypoint 

€£, fraction of available energy rate to be used for speed 
change in going to the Ith waypoint (discussed in 
appendix A) 

sines of the maximum and minimum inertial flightpath angle 

maximum and minimum true airspeed rates 

maximum and minimum values of energy rate to be used 

minimum landing speed with respect to the ship 

loiter speed (speed on the straight segment of the capture 
path) 

XAC,YAC,ZAC,HAC cartesian initial positions and course of the aircraft with 

respect to ship 


EPSLN(I) 

SGMAX,SGMIN 
VD GMAX , VD GMIN 
ETMAX , ETMIN 
VGLS 
VLN 
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VSH 


ship speed relative to Earth 


HSH ship heading relative to Earth 

Indicated airspeed is specified at each fixed-path waypoint. The speed-altitude 
profile along the synthesized horizontal path is obtained as follows: First, the 

equations of motion are integrated at constant altitude forward along the first cap- 
ture turn as indicated in figure 1. The speed changes at the specified maximum rate 
until (1) the loiter speed is achieved at A1 or (2) the turn is completed. In 
case (1) , the distance (Al - B) is stored and integration is continued to the end of 
the turn at B. In case (2), the integration continues until the loiter speed is 
achieved at A2 and the distance (A2 - D) is stored. 

Next, integration proceeds backward along the fixed path from the final waypoint 
to the capture waypoint with airspeed rate and flightpath angle specified, as dis- 
cussed later. If the speed and/or altitude specified for the next waypoint in the 
backward direction is reached before the waypoint, the distance traversed is stored 
and the speed and/or altitude change stops. If the waypoint is reached without 
achieving the waypoint speed and/or altitude, the speed and altitude actually 
achieved are stored and the target values are changed to those of the next waypoint 
in the backward direction. When the capture waypoint is reached, the target speed is 
set to the loiter speed and the target altitude is set to the aircraft altitude. 

When the backward integration reaches the end point of the forward integration 
(points A2 or Al) , the speeds and altitudes achieved at those points by the two 
integrations are compared* If they are equal, the synthesis was successful. If not, 
capture is not feasible along the horizontal path and a "no capture" message is indi- 
cated on the pilot f s display . 

Assuming a successful synthesis, a command table is generated which contains the 
following information for each state-change location: 

1. Mode flag (1 for speed and altitude changing, 2 for only speed changing, 

3 for altitude only changing, 4 for a constant-speed turn, 5 for straight, constant- 
speed flight , and 0 for no change) 

2. Index number of the next waypoint, assuming motion toward the final waypoint 

3. Correction factor used for reducing integration errors 

4. Time to go to the final waypoint 

5. Distance to go to the final waypoint 

6 . Flightpath angle 

7. Rate of change of airspeed 

8. Flap positions 

9. Bank angle 

10. Altitude 

11. Indicated airspeed 
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12. Lead distance for roll (compensates for finite roll rate) 

13. Lead distance for flightpath angle (compensates for finite flightpath angle 

rate) 

The information is used to initialize integration in real time forward along the next 
segment. Flight tests with NASA’s Augmentor Wing powered-lift airplane using a simi- 
lar strategy for the command table have shown that the aircraft can fly the trajectory 
defined by this table, reference 1. Modifications required to synthesize paths rela- 
tive to a moving ship are described in appendix C. 

The construction of the command table is performed with the help of the energy 
rate defined as the sum of aerodynamic flightpath angle and rate of change of true 
airspeed, E n = y a + V a /g. A detailed explanation of the technique is given in 
appendix A. The maximum and minimum energy rate at any flight condition characterizes 
the ability of the aircraft to change speed and altitude. Thus the algebraic sum of 
commanded values of Ya anc * V a /g must always be less than the maximum and greater 
than the minimum allowable energy rate. 

If a change in altitude or speed is called for, the system will attempt to make 
the change at the maximum values of sin Ya (which is found by converting SGMAX or 
SGMIN from inertial. to aerodynamic angle) or V a /g. If the aircraft is only . capable 
of an energy rate, E n , which is less than the desired value, the available E n is 
allocated according to 

V /g = eE 
a ° n 

sin y = (1 - e)E 
a n 

For simplicity, only three values of e are used, 0, 1, and 0.5. If e = 1 or 0, 
priority is assigned to (V a /g) or sin Y a » respectively. If j E n | exceeds the maxi- 
mum magnitude of the priority quantity, the remainder is assigned to the other. 

(For example, if e = l,JV a /g] < 0.05, and E n = 0.07, then sin y a = 0.02.) If 
e = 0.5, the ratios of V a /g and sin Y a to E n are maintained the same as those 
of the desired values . 

While the system will work with arbitrarily large magnitudes of desired sin Y a 
and V a /g, it is more efficient to set the operational limits to correspond as 
closely as practical to the aircraft’s capability. As the computer program is pres- 
ently implemented, small anomalies in the altitude profile or arbitrary; steps in Y a 
may occur at onset or at the end of a speed change where the. desired |E n | exceeds 
aircraft capability. The system provides smoothing of y a , V a /g, and altitude, but 
small "bumps” in Y a still appear in time histories. The values of e for each of the 
fixed waypoints is input with the waypoints, and all values for the capture path are 
set to a prespecified value. The value of VMAX for the capture waypoint is used to 
restrict the speed during the final capture turn and hence minimize the turn radius. 

REFERENCE FLIGHTPATHS 


Five reference flightpaths were used in this study. The flightpaths were based 
upon a TACAN approach chart for a prop or turboprop aircraft approaching a large, 
canted-deck carrier which is heading due north at 20 kt with its canted deck 
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centerlined at a heading of -11° • A 10-kt wind is blowing from -11°, or directly down 
the centerline of the canted deck. Initially, the aircraft is approaching the carrier 
along a course of 330° with respect to the canted deck (319° TACAN radial) with an 
airspeed of 140 kt (110 kt with respect to the ship) and at an altitude of 2000 ft 
(fig. 2). The distance to the touchdown point on the carrier deck is about 
12 DME n. mi. along the 319° TACAN radial. 


Path No. 1 

A horizontal profile with waypoints is shown in figure 2 for the first flight-* 
path. Path No. 1. This path uses a 1-n. mi. capture path to the first of the fixed 
waypoints that would be used by a pilot manually flying the path, which was taken 
from the Navy CV NATOPS MANUAL. The details of the waypoints for this path are given 
in table 1 wherein the position coordinates are in carrier-based coordinates in which 
the X-axis is positive along the magnetic north direction and the Y-axis is positive 
east. Each waypoint is specified by X-Y coordinates, an altitude, a turn radius, and 
a nominal and a maximum airspeed. This particular approach was chosen because the 
speed range of the Quiet Short-Haul Research Aircraft (QSRA) is more closely approxi- 
mated by a prop or turboprop aircraft. From a fuel and time point of view, this 
flightpath is more efficient than a manually flown flightplath because the ECG soft- 
ware system used to produce the reference path tries to minimize the fuel used along 
the capture path and between each fixed waypoint, whereas a pilot flying the same 
path could not normally be expected to be as efficient without the workload reaching 
an intolerable level. 


Path No. 2 

A profile with waypoints is shown in figure 3 for the second flightpath. Path 
No. 2. Details of the waypoints are shown in table 2. It differs from Path No. 1 in 
that waypoint 6 of Path No. 1 has been removed. The result is that the speed follow- 
ing the second turn is no longer forced to an intermediate 80-kt airspeed. Thus the 
ECG software system is free to maintain the higher 140-kt speed as long as possible 
before reducing the 65-kt speed required at 3 DME (slant range) n. mi. (waypoint 6 
of Path No. 2) . 


Path No. 3 

A horizontal profile with waypoints is shown in figure 4 for Path No. 3 which 
takes further advantage of the energy management capabilities of the ECG software 
system by greatly extending the length of the capture path and shortening the fixed 
path. Details of the waypoints are given in table 3. It used the same starting point 
as the previous two f lightpaths but was constrained to use only the same two final 
approach waypoints, that is, one waypoint at 3 DME n. mi. and 1200-ft altitude and the 
other at the touchdown point. 


Path No . 4 

A horizontal profile with waypoints is shown in figure 5 for Path No. 4 whose 
waypoint details are shown in table 4. It differs from Path No. 3 in that the way- 
point 3 DME n. mi. from touchdown was moved to 1 DME n. mi. from touchdown. This is 
a point on the final descent glide slope at which the airspeed is 65 kt and the 
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altitude is 919 ft. This flightpath saves additional fuel and time by further extend- 
ing the capture portion of the flightpath into the final approach where it is most 
productive. 


Path No, 5 

A horizontal profile with waypoints is shown in figure 6 for Path No. 5 whose 
waypoint details are given in table 5. This flightpath illustrates the ability of 
the ECG system to synthesize a capture path from any location. It illustrates best 
the full capabilities of the ECG software which, except for the constraints on the 
final 0.5 DME n. mi. of the landing approach, is free to exploit the maximum capabil- 
ities of the aircraft. The 0.5 n. mi. was chosen as a minimum distance which would 
allow an adequate amount of time for the aircraft to become stabilized on the final 
approach before touchdown. The initial conditions of the flightpath are the same as 
in the previous flightpaths except that the initial course is -90° and the X-component 
of position is the negative of the previous paths. Thus, the radial distance to the 
touchdown point in carrier-based coordinates is the same for all flightpaths. 

As for the previous paths, a knowledge of the aircraft’s location and velocity 
relative to the ship is assumed to have been provided to a sufficient accuracy by an 
on-board navigation system. This particular starting point was chosen because it 
provided the same radial distance from the ship as the other flightpaths and thus 
would allow comparison. 


RESULTS 


Described below are the results of the ECG synthesis software for each of the 
five reference flightpaths in response to a given set of initial conditions and fixed 
waypoints. In each case a sequence of actions is defined which, if implemented either 
manually or automatically, would fly the QSRA aircraft along the desired reference 
flightpath. 


Path No. 1 

Referring to figures 2 and 7, the QSRA aircraft proceeds along the 330° radial 
(relative to the canted deck) to waypoints 1 and 2 where a left turn (indicated by 
the roll attitude) is begun at a distance to go of about 86,000 ft. At the end of 
the turn, the distance along the intended path to the touchdown point is about 

79.000 ft or 13 n. mi. After level flight at an airspeed of 140 kt, the aircraft 
begins a constant-speed descent to an altitude of 1200 ft at a distance to go of about 

59.000 ft, and the pitch angle is reduced from about 2° to -0.5° during that period. 

A right, 60° turn is initiated to align the flightpath of the aircraft with the canted 
carrier-deck. As the turn progresses, the bank angle is reduced from 11° to 8° as a 
result of the ECG software trying to keep the turn radius constant while the air- 
craft’s speed with respect to the touchdown point decreases. After rolling out of the 
turn at a distance from touchdown of 6 DME n. mi., the aircraft resumes level flight 
at 1200-ft altitude along the centerline extension of the canted deck. The ECG soft- 
ware system then commands the deployment of the outboard flaps to a maximum of 59°, a 
condition required before the upper surface blowing (USB) flaps can be lowered. 

Later, at 28,000 ft to touchdown, the USB flaps are deployed to about 10° and the air- 
craft slows to 80 kt. When the aircraft has reached 19,000 ft to go, the USB flaps 
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are lowered an additional 26° to a total of 36° and the aircraft begins to slow to the 
commanded airspeed of 65 kt. As the altitude continues to decrease nearly linearly 
with distance to touchdown, the USB flaps are further extended to 43° at 7600 ft to 
go, where they remain until touchdown. The touchdown speed is approximately 35 kt 
with respect to the deck. The aerodynamic flightpath angle is -4.5°, which is the 
prespecified limit during the final approach to avoid an excessive sink rate at touch- 
down. The total fuel used during the descent is 1148 lb and the time required from 
the starting point is 747 sec. 

Effects of the various maneuvers on the amount of fuel used can be seen by exam- 
ining the plots of engine fan speed and fuel usage rate for each flightpath. 


Path No. 2 

The horizontal profile of this Path was shown in figure 3 and is similar to 
Path No. 1. As may be seen, by comparing figure 7 with figure 8, the two paths are 
the same until the aircraft has rolled out of the second turn at 6 DME n, mi. from 
touchdown and has resumed level flight at 1200-ft altitude. Figure 8 shows that the 
outboard flaps are deployed to 59° at a distance of 25,000 ft to go and the speed 
begins to decrease. This occurs about 7000 ft farther down the path than in figure 7. 
At about 23,000 ft to go, the USB flaps begin deployment and are out to 34° at about 
18,000 ft to go, at which time the aircraft’s airspeed has been reduced to 65 kt . 

The pitch angle is increased to 10.5° while the flightpath angle is being maintained 
at 0° since rollout from the last turn. At 7500 ft to go, the USB flaps are extended 
to 43° as the flightpath angle is established at -4.5° , and the pitch angle is main- 
tained at 6° where they are maintained until touchdown. The fuel used is 1080 lb and 
the elapsed time is 701 sec. 


Path No. 3 

This ECG flightpath has the same initial conditions as the previous two flight- 
paths, but the number of fixed waypoints is reduced to two. These are the two final 
approach waypoints from the set of six used in the previous path; that is, one way- 
point at 3 DME n. mi. and 1200 ft altitude and the other at the touchdown point. 

Figure 9 shows that upon leaving the starting point, the initial course is -40°, which 
is held for about 3000 ft. Then a turn is made to a course of -49° and is maintained 
until the aircraft nears the first of the final two waypoints. At 30,000 ft to go, 
the throttles are retarded slightly, and at 26,000 ft to go the aircraft begins to 
descend as the flightpath angle is established at -4.5°. At 24,000 ft to go, the 
outboard flaps are fully extended to 59°, and immediately thereafter extension of the 
USB flaps are started. Airspeed is about 105 kt . The USB flaps are extended to 44° 
at 22,000 ft to go while the airspeed has decreased to about 90 kt . The aircraft 
reaches an airspeed of 76 kt at 20,000 ft to go at which time the throttles are 
slowly advanced a small amount and the USB flaps are retracted to 31° to hold this 
speed while the aircraft banks 5° to make a right descending turn to a course of -11°, 
whereupon it levels out at 1200 ft altitude. During the course change, the speed is 
slowed to 65 kt . At 8000 ft to go, the flightpath angle is reduced rapidly from 0° 
to -4.5°, and the throttles are retarded slowly after a small step-reduction. The 
USB flaps are extended to 43°. The aircraft then continues in this configuration to 
touchdown; the elapsed time is 676 sec and the fuel used is 1070 lb. 
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Path No. 4 


The horizontal profile of Path No. 4 is the same as for Path No. 3 as may be 
seen by comparing figure 4 with figure 5. The only difference is that waypoint 1 in 
figure 5 has been moved 2 n. mi. closer to the touchdown point specified by way- 
point 2. This change gives the ECG software even greater flexibility to delay the 
final approach, and thus to exploit the aerodynamic capabilities of the powered-lift 
aircraft in a portion of the flightpath where the energy payoff is highest. Locating 
the first waypoint at 1 DME n. mi. from touchdown rather than at 3 DME n. mi. places 
it on the final glide slope approach at an altitude of 919 ft. The second waypoint 
is the touchdown point. 

Figure 10 shows that as the aircraft departs the initial position at an altitude 
of 2000 ft and at an airspeed of 140 kt, a small course correction is made at about 

83.000 ft to go, after which the aircraft reference path is undisturbed until about 

22.000 ft to go. At this point, the throttles are reduced and a slight pitch-down 
maneuver is executed, causing an altitude loss of about 300 ft. The airspeed then 
starts to decrease, followed shortly by start of the outboard flap deployment, which 
is complete by about 14,000 ft to go. At about 12,000 ft to go, pitch is commanded 
to about -5°, and the final descent begins with a flightpath angle which stabilizes 
at -4.5°. Also at about this time, the USB flaps are commanded to extend, stabiliz- 
ing at about 40°. As the USB flaps are extended, the pitch angle comes up to 6°. At 
about 7500 ft to go, a 5° bank angle is commanded to change the course from -44° to 
-11° to align the aircraft with the carrier's canted deck. A perturbation is seen in 
the USB flaps, throttles, and flightpath angle as the airspeed temporarily stabilizes 
at 75 kt during the final turn. At the end of the turn, the throttles are then 
advanced slightly, the USB flaps are set at 42°, and the airspeed is stabilized at 

65 kt as the aircraft continues down its descent path to touchdown on the carrier 
deck. Fuel used is 685 lb, which is a 385-lb reduction over the previous path. The 
elapsed time is 529 sec. 

The fuel use reduction for Path No. 4 relative to Path No. 3 can be explained 
with the aid of figure 11. The synthesized data are plotted against time to compare 
the fuel use and time spent at each flight regime for the two approaches. Path No. 3 
requires a 676-sec flight time whereas Path No. 4 requires only 529 sec, a savings of 
147 sec. Figure 11 shows that the fuel rate is the same (about 1.3 lb/sec) for both 
paths as long as the airspeed is maintained at 140 kt. This is held for 281 sec in 
Path No, 3 and for 335 sec, or 54 sec longer, in Path No. 4. At these times, the 
engine speeds are reduced in both cases to begin decelerating with a corresponding 
drop in the fuel-flow rate. As the airspeed slows below approximately 100 kt, the 
outboard flaps are extended and the throttles reduced again. When the outboard flaps 
are fully extended, the USB flaps begin extension, causing the airplane to enter a 
powered-lift mode during which engine thrust is being used increasingly for lift, 
necessitating an increase in engine speed and thereby causing the fuel-flow rate to 
increase. Net fuel usage, however, is less than before on both paths during the 
deceleration and descent from the 2000- ft altitude phase. Path No. 3 then has a seg- 
ment of about 165-sec duration at a constant altitude of 1200 ft which is not in 
Path No. 4. During this period, the fuel rate for Path No. 3 is at its highest 
level, between 1.8 and 2.3 lb/sec. The increase during this interval is caused by 
more engine thrust being required to hold both the altitude and speed constant as the 
USB flaps are being slowly retracted then extended again for the final approach. 

Thus, on Path No. 4, the final glide slope was captured while the descent from 2000 ft 
was still in progress, but Path No. 3 had a constant altitude and airspeed segment, 
which had a high fuel rate. The paths and the fuel rate are the same during the final 
100 sec after the common glide slope was reached. The savings in fuel for Path No. 4 

9 


j 



can, therefore, be traced to a combination of (1) a shorter total flight time; 

(2) flying 54 sec longer at 140 kt , and eliminating 165 sec spent holding altitude 
and airspeed in aircraft configurations in which significant engine thrust is 
required to enhance lift; and (3) a continuous transition from descent to final 
glide slope. 


Path No. 5 

The starting point of this path has a radial distance from the carrier which is 
the same as the previous paths, but the distance along the path is about 10,000 ft 
greater than Path No. 4. It is also different from the previous paths in that the 
initial course is -90° and the initial X-position is the negative of the previous 
paths. A horizontal profile of the path is shown in figure 6. This figure shows that 
the first waypoint is on the final approach at 0.5 n. mi. from touchdown. This gives 
the ECG software the maximum flexibility to save fuel and time while giving the air- 
craft adequate time to become stabilized on the final approach. 

Figure 12 shows that as the aircraft departs the starting point at an altitude 
of 2000 ft and at an airspeed of 140 kt, a 10° left bank begins at about 94,000 ft to 
go. The bank angle is increased as the aircraft's speed relative to the landing point 
increases as a result of the ECG software trying to keep the turn radius constant. 

The aircraft rolls out of the turn on a course of -149° and continues at the initial 
speed and altitude until about 32,000 ft to go. The aircraft is pitched down slightly 
as the throttles are retarded slightly, causing the aircraft to reduce its airspeed 
and begin a descent. At about 23,000 ft to go, the outboard flaps begin to be 
deployed to 59°. When the aircraft has reached about 20,000 ft to go, the outboard 
flaps are fully extended. Then, the USB flaps begin to extend at the maximum rate 
while the pitch angle and angle of attack are increased positively until about 

10,000 ft to go. The final turn to align the aircraft with the canted deck begins and 

the flightpath angle is established at its final value of -4.5° for the remainder of 
the approach, but the airspeed is high at about 75 kt. The USB flaps are lowered from 

30° to 42° at about 3000 ft to go and the throttles are increased slightly as the 

airspeed is reduced to final speed of 65 kt . The fuel consumed is about 552 lb and 
the time is 442 sec. 


CONCLUDING REMARKS 


A summary of some of the more important results of the flightpath synthesis for 
the five approaches described in the preceding section are shown in table 6. Dis- 
tance flown, flight time, and fuel used are presented for each synthesized flight. 

Note that, as previously described, the first four approaches started at the same 
location, course, and airspeed. Also, all four approaches had the same course, air- 
speed, horizontal profile, and glide path from 1 DME n. mi. to touchdown. The fifth 
path started at a different location, but flew the same horizontal profile, glide 
path, and airspeed as the other four from 0.5 n. mi. to touchdown. The differences 
between the five approaches illustrate the unique capabilities of the ECG system. 

As shown in table 6, the baseline Path No. 1 approach flew 95,062 ft to touchdown 
in 747 sec and used 1148 lb of fuel. Although Path No. 2 flew the same predefined 
horizontal and vertical profiles, allowing the ECG system to fly a higher speed for a 
longer period shortened the flight time by 46 see (a 6.2% reduction) and saved 68 lb 
of fuel (5.9%) less). When the ECG system was also allowed to exploit its horizontal 
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and vertical flightpath synthesis capabilities to capture a waypoint nearer the 
touchdown point, the Path No. 3 flight was completed in a 71-sec shorter flight time 
(9.5%) and used 78 lb (6.8%) less fuel than the baseline approach while still provid- 
ing the same stabilized final descent as the other two approaches from 3 DME n. mi. 
Even more fuel can be saved by capturing a waypoint which is much closer to touchdown. 
This is seen in Path No. 4, in which the captured waypoint was on the final glide 
slope only 1 n. mi. (slant range) from touchdown. This procedure saved 218 sec 
(29.2%) of flight time and 463 lb of fuel (40.3%) relative to Path No. 1. These 
savings were produced by (1) shortening the total flight time and (2) keeping air- 
speed high as long as possible to minimize time spent flying in a regime in which 
more thrust is required for lift to aid the aerodynamic lift, a less fuel efficient 
condition. Significant fuel savings were also shown when time spent flying at con- 
stant altitude at low airspeeds was minimized or eliminated. 

The fuel savings for Paths No. 2, 3, and 4 indicate the potential of the ECG 
system when there is some freedom to select the horizontal, altitude, and speed pro- 
files. More fuel could be saved in situations in which the system’s path generation 
and flight time prediction capabilities could be used to coordinate direct landings 
by aircraft without spending time in holding patterns, as indicated by the results for 
Path No. 5. As the ECG system is allowed more freedom, fuel consumption decreases. 

Use of the total-path-synthesis procedure used in Paths 3, 4, and 5 may depend upon 
policy changes to the extent of allowing more flexibility in the selection of the 
approach path between the holding pattern and the final descent on the glidepath. 

Even if that flexibility is not allowed for operational, communications, or security 
reasons, this study, although limited in scope, shows that an ECG system can save 
flight time and fuel by more efficiently managing the aircraft’s capabilities. 

Defining operationally acceptable limits on such parameters as minimum allowable 
flight time on a stabilized glide slope to touchdown and the allowable horizontal, 
vertical, and airspeed profiles would require experimental tests of the ECG concept. 
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APPENDIX A 


ENERGY-RATE EQUATIONS AND SIMPLIFIED AIRCRAFT MODEL 


The basic objective of the energy -management guidance system is to synthesize a 
flyable trajectory from the current aircraft states to some desired terminal states. 
The term "flyable" implies compliance with certain constraints. 

The first set of constraints (which will be referred to as operational con- 
straints) are externally imposed and are specified by input constants. Some opera- 
tional constraints, including limits on inertial flightpath angle, airspeed rate, 
normal acceleration, bank angle, and roll rate, are dictated by passenger comfort or 
pilot preference. Others, determined by Air Traffic Control regulations or by 
terrain, take the form of specified horizontal paths together with speed and altitude 
restrictions. 

Next are geometric constraints, which require that maneuvers meeting the opera- 
tional constraints can actually be made in the space available. For example, given a 
maximum flightpath angle |Y m | 9 a maximum rate of change of |y| (hence a maximum 
normal acceleration), and a desired change in altitude Ahp, is it possible for |y| 
to change from its present value to |Y m | and then to its final value while Ahp is 
being accomplished? Where should pitch-over to the final value of y begin? These 
problems are solved by the system through closed-form equations which will be 
described in appendix B. A similar problem for the lateral case is handled by intro- 
ducing step changes in the reference bank angle. 

The operational and geometrical constraints are only indirectly affected by air- 
craft performance and do not guarantee that the aircraft can perform all maneuvers 
which satisfy them. The final set of constraints (which will be referred to as air- 
craft constraints) restrict the synthesized flightpath to flight conditions and con- 
trol settings which the aircraft can achieve while meeting safety standards and 
reserving maneuver capability for control. Determination of the aircraft constraints 
requires a fast-time solution of the aircraft’s equations of motion. The algorithms 
used for this solution depend upon energy-rate methods and are described in this 
appendix. 


Energy-Rate Equations 

The aircraft energy is defined as the sum of the kinetic and potential energies 

E = mgh + mV 2 (Al) 

z a 


where h is the altitude. 

Differentiating equation (Al) with respect to time and dividing by mg gives 

| = h + — V (A2) 

W g a 
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Substituting V a sin y a for h and dividing by V a gives the normalized 
energy rate E n as 


E 

n 


E 

WV 

a 


V 

sin y 4- — 
a g 


(A3) 


For convenience E n will be referred to simply as the energy rate* 

It can be seen from equation (A3) that E n represents the total capability of 
the aircraft to change speed and altitude. Given a desired flightpath angle and air- 
speed rate, it is necessary to find the control settings which will produce the 
desired E n in steady state. 


Equations of Motion 

The trajectory-synthesis algorithm requires a simple, compact method for inte- 
grating the longitudinal equations of motion while meeting all of the imposed con- 
straints (operational, geometric, and aircraft). This is accomplished through the use 
of the following point-mass equations of motion: 

mV = T cos a - D - RAMDRAG - mg sin y - mV cos y (A4) 

a a wat 'a 

mV y^ = (T sin a + L)cos $ - mg cos y - mV ^ sin y (A5) 

a a a wat a 


where 


T 

D 

RAMDRAG 

L 

mg 

Y a 

a 

V _ 
wat 


total gross thrust 
aerodynamic drag 

drag from ram air through engine 
aerodynamic lift 
aircraft weight (W) 
aerodynamic flightpath angle 
bank angle 
angle of attack 

along-track component of rate of change of the mean wind 


It is assumed that the mean wind is a function only of altitude and has no vertical 
component so that 


av 


V 


wat 


av 


wat 


ah 


h - 


wat 


ah 


V sin y 
a a 


(A6) 
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where V wat is the along-track component of the wind velocity and h is the 
alt itude . 


The lateral equations of motion per se need not be integrated, since the effects 
of bank angle are accounted for by using a fictitious aircraft weight equal to the 
true weight divided by cos cf> . 


The terms in equations (A4) and (A5) are converted to normalized coefficient 
form by dividing by dynamic pressure, q, and wing area, S w , to give 


~ ~ ~ (Cp. + C ) - sin y - V cos y 

qS w g D RD 'a wat 'a 


(A7) 


W V a • p 

^ Y = C T cos 

qS g a L 
’ w & 


W , , wat . N 

(cos y H — - sin y ) 

qS w a g V 


(A 8 ) 


where 


C = -(T cos a - D)/qS 


w 


and 


= RAMDRAG/qS w 


C T = (L + T sin a)/qS 

L W 


Multiplying equation (A3) by W/qS w and substituting from equation (A7) gives 

V 

W v _ fn r- \ W wa t 

qS w E n " (C D C RD qS w g COS 7 a 


(A9) 


Equation (A 8 ) requires 


C = — 


W 


L qS cos 6 
w 


V 

. wat 

Y + cos y -4 sxn y 

a a g a 


) 


or, since for trim y = 0 , 

a 


— — r 


W 


L qS cos 
n w 


(■ 


V ^ 

, wat 

cos y + sm y 

a g a 


) 


(A10) 


Solving equation (A9) for -(Cp + Crx>) gives 


V 

tn , n \ W ^ , W Wat 

■(G n + C ) = E + =5 cos y 

D RD qS n qS g a 


(All) 
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(A12) 


The coefficients Cp, Crd, and Cl are functions only of <%, Op, and flaps (<5 fref) 


Energy-Rate Tables 

The values of Ex, the corresponding wing angle of attack, a w , the thrust coef- 
ficient, Cx» and the flap parameter, <SfrEF> were calculated and stored in what are 
referred to as energy-rate tables. The values of from these tables are plotted 

in figure A1 as a function of lift coefficient. Cl, for a set of engine parameters 
and weights defined in table Al. The indices in table A1 correspond to the curves in 
figure Al in sequence from the uppermost curve. Note that the pressure ratios 1.0, 
0.688, and 0.459 correspond to pressure altitudes of sea level, 10,000 ft, and 
20,000 ft, respectively. Altitude is accounted for by linear interpolation on pres- 
sure ratio. The nominal aircraft weights listed in table Al were used to establish 
airspeed, which, in turn, determines CxMAX* and the placard values of flaps. 

For the minimum aircraft weight (40,000 lb), and the placard flaps are the 

maximum value for any given value of for maximum weight (56,000 lb), CpjjiN is 

the minimum value for each C L . The flap placards for the minimum aircraft weight 
were used for generating all the energy-rate table data. Curve 1 in figure Al repre- 
sents the maximum E T atainable for any weight and altitude while curve 7 represents 
the minimum. Thus, those two curves represent the limits to the flight envelope over 
the range of Cl« The effects of spoilers and various other complications which 
affect the algorithms used in the program are discussed in the next section. 

During synthesis, desired values of sin y a and V a /g which satisfy the opera- 
tional and geometric constraints are obtained and used to calculate the desired values 
of E n . The desired value of E n together with its operational maximum and minimum 
are used in equation (A12) to find the desired, maximum, and minimuih values of Ex« 
Then linear interpolation on the current value of Cl is used to obtain a one- 
dimensional array of E^. This is equivalent to finding Ep at the intersection of 
a vertical line on figure Al at the current value of Cl with Cp from the different 
thrust contours. Arrays of a w , Cx» and flap settings corresponding to this value of 
Cl are also found. The constraints of thrust and flaps are then applied to obtain 
the maximum and minimum values of Ep achievable. These are stored in the first and 
last elements of the Ep array, respectively, and the corresponding values of o^, 

Cp, and flaps are adjusted accordingly. Note that where the contours merge, the 
number of elements in the one-dimensional arrays is reduced. Finally, the maximum and 
minimum values of Ep are multiplied by a factor of a (currently 0.9) to reserve the 
fraction (1 - o) of the available Ep for control. If the desired value of Ep com- 
puted from equation (A12) exceeds the limits on tabular values of Ep just discussed. 
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it is set equal to the tabular limits and equation (A12) is inverted to find the 
resulting limit on E n . 

After the desired value for E n has been determined to be within operational 
and geometric limits, sin y a and V a /g are recomputed using an algorithm controlled 
by the pilot-input parameter e to determine the amount of E n allocated to each. 
Because E n is a function of sin y a , iteration is used to find the final values of 
E n , y a , V a /g, a w , thrust coefficient, and flaps. 

For a given set of a w , C^, and 6 fref> the same value of E n will result, inde- 
pendent of speed or altitude. However, the maximum value of <$FREF which can be 
used (the placard value) is a function of airspeed. Since, for trim, it is necessary 
that ~ W/qS w , the maximum value of 6 frBF is a function of aircraft weight. The 
values of E^ were calculated using maximum and minimum aircraft weights as discussed 
earlier. The placard flaps for the minimum weight were used for these tables and are 
the maximum that will be encountered. For heavier weights the maximum allowable 
^FREF will be reduced, resulting in turn in the reduction of the maximum magnitude 
of negative energy rate achievable. Interpolation on ^FREF is used to adjust E»£ 
to a value achievable for the actual aircraft weight. 

Calculation of Energy-Rate Tables - If we assume that we have a desired flight- 
path angle and airspeed rate, we wish to determine whether the desired E n can be 
achieved at the current flight conditions, and if so, to determine the corresponding 
required control settings for flaps, angle of attack, and thrust. Sets of tabular 
data referred to as ’’energy-rate tables” have been generated for this purpose. 

The first step in generating the energy-rate tables was to compute sets of static 
trim conditions for the longitudinal mode using the CDC 7600 simulation software of 
the QSRA. The simulation uses tabular data from wind tunnel experiments described in 
reference 5. The simulation was used to compute sets of tables of C i and E>p versus 
a w for constant values of Dj and flaps. Here C<£ is defined as the total gross 
thrust for four engines divided by qS w . This operation produced sets of tables, 
examples of which are illustrated in graphical form in figure A2. Tables representing 
the various combinations of outboard flaps (6y) , USB flaps (^uSb) 5 an< ^ % sufficient 
to cover the aircraft flight envelope were generated. Values of C«£ and a w used for 
the original wind tunnel data were also used in this calculation. 

Note that the flight control system is designed so that <$USB constrained to 
0° unless = 59°. Hence, and may be regarded as a single control <$FREF* 

Also note that it is assumed that the spoilers are biased from 0° to 10° proportion- 
ally to the value of 6usB between 0° and 30°. As will be discussed later, the data 
indicate that it may be necessary to increase the spoiler bias in order to obtain 
sufficiently negative values of E n . 

For the synthesis process, it is necessary that be an independent variable, 

and it was necessary to "invert” the trim data using linear interpolation. Also, 
since there is one more control than necessary to trim the aircraft, some criterion 
must be used to eliminate the redundancy, even in the case of manual operation. 
Furthermore, the elimination process must be implemented on the computer if a refer- 
ence flightpath is to be synthesized. The criterion chosen for the initial applica- 
tion to the QSRA is the minimization of fuel. A set of energy-rate tables were gener- 
ated using the minimum-fuel criterion by ^ finding, for specified combinations of 
and C T , the value of SpREF for which E-j- is maximum. The minimum feasible value of 
E n was also found, assuming that = 0 and with the constraint that 

0 £ a w < 15°. Placard values of 6 freF and C TMAX were computed for an aircraft weight 


4 


16 



of 40,000 lb. Note that when Cl is calculated for the interpolation of the Cl 
table described in the section on speed-altitude synthesis, the value of y a in 
equation (A10) was set to zero. This was done to reduce computation time, but if more 
precision is required, iteration would be necessary. 

The results of these computations in figure A3 show as a function of 

for different values of Cy. The upper dashed line is the contour of the maximum 
allowable value of Cj resulting from the constraint that the corrected engine core 
speed, N 2 //§7, be less than or equal to 96%. This computation was made assuming 
6 am B = 1 (sea level on a standard day). 6 ^mb only rarely exceeds unity, but 
decreases exponentially with increasing altitude. CtMAX ma y be further reduced as 
the result of the constraint that the engine exhaust gas temperature not exceed 920° C. 

The lower dashed curve is for the minimum value of E«j obtainable with the 
spoiler values specified earlier and no minimum constraint on C^. For values of Cl 
greater than that (Cl = 2.4) at the intersection of the lower dashed line and the 
Cx = 0 contour, the value of E n is constrained by = 15°. At lower values of 
Cl, the constraint is the result of the flap placards. The two dashed curves, 
together with the vertical lines at the maximum and minimum values of Cl> represent 
an outer limit to the operational envelope of the aircraft. The maximum value of 
Cl = 7 results in minimum equivalent airspeeds of 53 kt and 63 kt for aircraft 
weights of 40,000 lb and 56,000 lb, respectively. The minimum value of Cl corre- 
sponds to the maximum equivalent airspeed of 160 kt and an aircraft weight of 
40,000 lb. The minimum value of Cl will increase with aircraft weight. 

Figures A3(b) and A3(e) show the values of a w and 6pREF» respectively, corre- 
sponding to the values of 'Ey in figure A3 (a). The rather erratic nature of these 
curves resulted in the addition of a number of intermediate values of <$freF to t be 
trim data before computing the energy-rate tables finally used, which are shown in 
figure Al. 

Figure A4(a) shows curves of Ej versus C L obtained by interpolation of data of 
the type used for figure A3. The upper three solid curves are for C^MAX at va l ues 
of 6^ B corresponding to pressure altitudes of sea level, 10,000 ft, and 20,000 ft. 
The values of E^ shown are the maximum obtainable with those values of C^. The 
lower two solid curves are for the maximum E T obtainable using the minimum allowable 
value of Cp (N 2 //e^ - 78%) at sea level and 10,000 ft. The two dashed curves are the 
minimum values of E n obtainable with the minimum values of Cp. The dotted curve 
is the minimum Ep obtainable with CpMiN set to 0. Note that in contrast to 
figure Al discussed in the beginning of this appendix, minimum aircraft weight 
(40,000 lb) was used for all of these curves. This choice simplifies the following 
discussion. However, the data actually used in the program are represented by 
figure Al. At the higher values of Cl, all of the five lower curves merge into a 
single curve determined by the constraints on a w and 6 ]?ref- 

Figures A4(b), A4(c), and A4(d) show the corresponding values of a w , 6p^£p» anc * 
Cp, respectively. In all of the maximum Ep cases (solid curves), rises linearly 

with Cl until 6 rreF changes from 0% to 30%, At this point, a w decreases to a 
lower value and then increases linearly until it reaches 15° (the maximum allowable 
value). Some of the curves have a break point and decrease in slope near a w = 15°. 
This change in slope could be eliminated by finding the exact value of Cl for which 
a w attains 15°. Also, 6p R g F tends to remain at 30% until a w reaches 15°. Similar 
relationships exist for the minimum Ep curves (dashed lines) , but are complicated by 
the flap placards and by the nature of the trim data at high values of 6ugg. 
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Figure A4(d) shows the thrust coefficients which correspond to the data for 
figure A3(a,b,c), There are only five curves since both the minimum and maximum 
energy-rate data are plotted for C^mIN at sea level and 10,000 ft. At lower values 
of C^, the values of are linear functions of C^. The changes in slope at the 

higher values of C;l are caused when the limits on Sjggp and a w are encountered. 

It can be seen from figure A4(a) that the descent /deceleration capability of the 
QSRA is rather limited. For > 2.0, the minimum value of Ex at sea level has an 

average value of about -0.15. For higher weights, the minimum value of Ej may be 
larger in magnitude depending on the flap placards. However, it is desirable to 
increase the descent /deceleration capability, particularly at the lower weights and 
smaller value of Cl* One possible method of improving the descent /deceleration 
capability was investigated, and the results are presented in figure A5. 

Figures A5(a) and A5(b) show the minimum Ej curves for G^iN (<$AMB = 1) an< ^ 

CtMIN (^AMB = 0.687) from figure A4 (a) as solid lines. The dashed curves are the 
result of increasing the spoiler bias by 8° from the nominal for all values of Cl- 
(The nominal spoiler bias increases linearly with USB flaps at the rate of one-third 
of a degree per degree of USB to a maximum spoiler bias of 10°. ) This change pro- 
duces a significant effect on the minimum E^ for C^ < 2.0. The effectiveness of 
the spoilers in producing more negative values of E^ is reduced with increasing 
6 usb an< i to a l esser extent with increasing C^. 

It is apparent from the results in figure A5 that this increase in spoiler 
deflection is not effective in producing more negative values of E n much beyond 
Cp, = 2.0 (i.e. , for equivalent airspeeds less than about 100 kt) . (It has been sug- 
gested by NASA personnel that even larger increases in spoiler deflection are feasible 
and may produce desirable results. Further study of the problem is needed.) 


Engine Speed Constraints on Thrust Coefficient 


The thrust coefficient, is constrained to a maximum value CxMAX 5 w hich 

requires that the corrected engine core speed, N 2 //0^, not exceed 96% rpm while a 
lower limit, C^iN > is established by the requirement that N 2 //oJ > 78% rpm. The 
values of the resulting limits on Tq/ 6^ s , where T q is the gross thrust for four 
engines, are given in table A2 . 


Mach number is defined as the ratio of true airspeed, to the speed of sound, 

a. By definition of VgQ, the equivalent airspeed 



(A13) 


where p is the actual density and p 0 is the sea level density of the standard 
atmosphere (ref . 6), It can be shown p/p Q = ^ AMB / ^AMB ^h^re Q^MB is the ratio of 
the ambient absolute temperature to the standard value at sea level. 


It can also be shown that 


a/a 

o 



(A14) 


where a Q is the speed of sound at sea level in the standard atmosphere. Substitut- 
ing equations (All) and (A14) into the definition of Mach number, M, gives 
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(A15) 



1 



Linear interpolation on Mach number is used to find TqmAx/^AMB anc ^ ^GMIN^AMB’ 
which are then multiplied by <$AMb/5S w to give ^XMAX an ^ ^TMIN • 


Calculation of Limit on Cx Due to Exhaust Gas 
Temperature Limit of 920° C 


The data used for these calculations were obtained from CDC 7600 simulation 
tables TMGTl (HPBLC =0) , TMBT2 (HPBLC=5) , and TMGT3 (HPBLC-10) , where high-pressure 
boundary-layer control (HPBLC) is the percentage of cold thrust air bled from the 
engines to be blown over the leading edge and the ailerons for boundary layer control. 
These are tables of MGT / 0 2 * 8 8 1 versus and M, where MGT is the measured 

(exhaust) gas temperature, M is the Mach number, Nj is the fan speed (in % rpm) , and 
0 2 is defined by 0 2 = (1.0 + 0 . 2M 2 )T^xb/T 0 . Here T^mb is the ambient absolute 
temperature and T 0 is the standard value of T&ixg at sea level, or we can define 
t AMB ® T o + AT. By using these tables it is possible for a given value of 0 2 to 
find the value of N 1 //0^ for which MGT = 920° C. When this was done for values of 
AT in the range -30° C < AT < 30° C and for 0 < M < 0.3, it was found that Mach 
number has a negligible effect on the value of ^i/* / 0^"- That is, if we define Nxmx 
as the value of for which MGT = 920° C, then 


N 

TMX 


N + kAT 
o 


where N Q and k are functions of HPBLC. A conservative approximation is given by 

N_ mv = 94.0 - 0.78 (HPBLC) - AT[0.1667 4- 0 .00333 (HPBLC) ] (A16) 

1 MX. 

This result is plotted in figure A6 for the three values of HPBLC for which data are 
available. Note that Nxmx strongly dependent on HPBLC. The Boeing math model 
report (ref. 3) states that HPBLC varies from 10% at STOL idle to 0% at high thrust 
values . 

The value of HPBLC in the CDC 7600 simulation is the result of a set of calcula- 
tions using the actual core speed N 2 (not N 2 //0^ ) • Given a value of N 2 , the value 
of HPBLC is adjusted so that the temperature constraint is not violated unless 
N 3.//07 is greater than the maximum for HPBLC = 0, in which case HPBLC is set to 0. 
Note that HPBLC is not necessarily adjusted to its allowable maximum for the given 
N 1 //0^", but never exceeds it. Therefore, the maximum N x / / q^ for MGT < 920° C 
can be taken as the value for HPBLC = 0, that is, 

N x //0^ = 94.0 - AT/6.0 (A17) 

The corresponding value of Tq/6 amb can f° un d in table A3 and multiplied by 
<$amb/5 s w to £ et t ^ ie max i mum value of thrust coefficient , Cxx* which satisfies the 
temperature constraint. Then Cx^AX as set to the minimum of its original value 
(for which N 2 //0^ < 96%) and C^. 
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APPENDIX B 


DERIVATION AND USE OF CLOSED-FORM EQUATIONS FOR 
SPEED-ALTITUDE PROFILE 


In order to determine the aircraft’s capability to achieve a desired flightpath 
angle during a predetermined altitude change, a set of closed-form equations has been 
developed which calculate the altitude changes during changes in y and/or V. When 
y changes, its maximum magnitude is subject to the constraint that ]y| - k^/V^, 
where Vj is the true airspeed and k x is an input normal acceleration. Note that 
if speed and altitude change in opposite directions, the total change in energy is 
minimized by requiring that the change in energy caused by speed change and that 
caused by altitude change have a fixed ratio which will cause them to be completed 
simultaneously. Four equations are used in the calculation of these altitude changes: 
one for V and h changing in the same direction, one for V constant while y 
changes, one for y constant and V changing, and one for V and h changing in 
opposite directions. These equations are applied in sequence to the intervals of the 
trajectory between successive pairs of waypoints during which the altitude changes. 

Case 1: V and h change: same sign 


We assume 


V = V 0 + k 2 t 

Y = kj./v 

h = V sin y 


From equations (Bl) and (B2) 


k i 

v 0 + k 2 t 



“ ^1 



dt 

V 0 + k 2 t 


Integrating equation (B5) gives 


Y 



ln(y n + k,t) 


or 



(Bl) 

(B2) 

(B3) 

(B4) 

(B5) 


(B6) 


(B7) 


Let u = y - Y 0 * We have 


JS 
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kl / k 2 \ 

u 'k; lB V 1 + n V 

Inverting equation (B8) 

T [e (k2u/kl) - 1] 

k 2 

and differentiating from equation (B9) 

dt - ^ e (k 2“ /kl) du 

k i 

Substituting equation (B9) into equation (Bl) gives 

V - V„ g( k 2 u / k l) 

and substituting equations (BIO) and (B 1 1 ) into equation (B3) gives 

dh = e ( 2k 2 u / k i) sin(y 0 + u)du 

k i 

Let Vo/k x = b and 2k 2 /k 1 = a for convenience. Then 

dh = b e au sin(v 0 + u)du 
Using double angle formulas with (B13) gives 


dh = b e au [sin y 0 cos u + cos y Q sin u]du 


so that upon integration. 


h - h Q = b ^sin y 0 [e au (a cos u + sin u) ] ^ a f 

au 1 Vl u 

+ cos Y 0 [e (a sin u + cos u) ] \ + a a j | 0 


or 


(B8) 


(B9) 


(BIO) 


(Bll) 


(B12) 


(B13) 


Ah = h - h n = 


b {[a sin(u + y 0 ) ~ cos(u + Y 0 )]e au - [a sin y 0 - cos y q ]> 


(B14) 


o 1 + a 2 

Substituting y - y 0 for u in equation (B14) gives 

= l — ~ — 2 [e a ^ Y_Y °- ) (a sin y - cos y) “ (a sin y 0 - cos y Q )] (B15) 

jl t a 


Ah 
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Factoring e 


from equation (B15) gives 


-ay 0 


, _aY 0 

Ah = — — [e a ^(a sin y - cos y) - e aY °(a sin y Q - cos Y 0 )] (B16) 

1 + a 2 

where Ah is the altitude change while changing y and V. 

Case 2: V constant, h changing 


In this case we desire to specify Vgq instead of Vj because it is equivalent 
airspeed which will be held constant. A good approximation for the altitudes under 
consideration is 


V T . V EQ (1 + Kh) 

where K = 0.16xl0~ if ft"" 1 and h is in feet. 

From equation (B2) 


Y " V I ' V EQ (1 + KW 


h = V„ sin y = V__(l + Kh)sin y 
1 JiLJ 


From equations (B 1 7 ) and (B18) 


dh ^ fv E q(i 4 


dy 


ki 


sxn y 


so that 


or 


ki dh 


v“ Q (i + Kh ) 2 


= sin y dy 


V. 


EQ •'ho 

Integrating equation (B19) gives 

k 


i. r — ( 

„ 1 (1 + Kh) 2 I 


sin y dy 


r i — _ 1 1 = 

2 [l + Kh 0 1 + KhJ 
EQ 


cos y 0 - cos y 


Solving equation (B20) for h results in 


(B17) 

(B18) 


(B19) 


(B20) 
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(B21) 


h= i 


k^d + Kh 0 ) 


-K(l + Kh 0 )V EQ (cos y o “ cos y) + k 




From which it follows that 


Ah = h - h Q = 


(1 + Kh 0 ) 2 Vg Q (cos Y Q - cos y) 
k x - K(1 + Kh 0 )V EQ (cos Y o " cos y) 


Case 3: y constant, V changing 

From equations (Bl) and (B3 ) , 

h = (V 0 + kg t) 'sin y 0 

V -k 2 

dh = •£- sin y 0 dV 


Integrating equation (B23) gives 

Ah = h - h 0 



(B22) 


(B23) 


Case 4; V and h changing in opposite directions 

In order to minimize the total change in energy (and thus minimize |E n |), we 
desire the altitude and speed changes to be accomplished simultaneously for this case. 
We define the aircraft energy E as 


E = h + (B24) 

2g 

Therefore, 

AE, = h - h n - Ah (B25a) 

h 0 

V 2 - V* 

AE = 7T — - (B25b) 

v 2g 

where AE^ is the change in energy from altitude change, and AE V is the change in 
energy from velocity change. Differentiating equations (B25a) and (B25b) , 

(B26a) 


(B26b) 


AE h -h 


• w 

AE = — 
v g 
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Define 


AE AE 
k* = v 


AE. 


h AE, 


Then to satisfy the minimum total energy change condition, we want 

VV 

gh 


w -k* 


or, substituting V sin y for h, 

V ® k*g sin y 

or, defining K 3 = gk*, 

V = K 3 sin y 

We want to restrict y to y = ki/V (from eq. (32)), so that 

dv V k 3 Sin Y 


from which 


dy y ki/V 


dV k 3 

-y- = k" sin y d V 


and, integrating equation (B32) 


f 

w n 


dV 3 

T “ “ k^ cos y 


V 

in — - — (cos y q - cos y) 
v o 


V * V 0 e 


(k 3 /ki) (cosy 0 -cosy) 


dh _ V sin y 
dV k 3 sin y k 3 


and, integrating equation (B34) 


Ah = h - h n 




i 


(B27) 

(B28) 

(B29) 

(B30) 

(B31) 

(B32) 

(B33) 

(B34) 
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Computer Implementation of Closed-Form Equations 

The closed-form equations developed in this section are implemented in the com- 
puter program subroutines INTEG and VGMALT described in appendix C. There are two 
sections in INTEG which solve the equations for: (1) Case 4, where the speed and 

altitude changes called for are of opposite signs, and (2) Cases 1, 2, and 3. 

Case 4 is solved in a simple iterative loop. Since the velocity-altitude rela- 
tionship is known in this case throughout the interval, it is possible to determine 
the altitude changes during the required changes in y at both ends of the interval 
in closed form. 

After this is done it is determined whether the sum of these altitude changes 
exceeds the total altitude change desired in the trajectory; if so, y c is limited to 
its old value multiplied by the ratio of the total altitude change desired over the 
sum of the individual altitude changes, and the computation is iterated. 

It is possible that the altitude-velocity relationship may not be solved in one 
calculation at the end of the trajectory for Cases 1, 2, and 3, since finding the 
velocity at the beginning of the final pitch maneuver may require an iteration pro- 
cess. Therefore, the process used for these cases is somewhat different from that 
for Case 4. 

At the beginning of the closed-form process for Cases 1, 2, and 3 it is deter- 
mined whether or not there is to be a velocity change during the interval. If not, 
program execution skips directly to the constant equivalent airspeed portion of the 
logic. If there is a velocity change, the subroutine VGMALT is called. 

The first step in VGMALT is to solve equation (B9) for the time, ty, required to 
achieve y c and then to solve equation (Bl) for the time, t v , required to achieve 
the desired speed change. If t v < ty* equation (8) (with x - t v ) is used to com- 
pute u, the change in y, and if t v > ty, u is set to y c - y Q . Equation (B16) is 
then used to compute the corresponding altitude change. Ah, which is stored in DELHG 
if the y change is completed first and in DELHV if the speed changed is completed 
first . 

If the speed change is completed first, program execution continues to the con- 
stant equivalent airspeed portion of the closed-form algorithm, where equation (B22) 
is used to compute and store in DELHG the remaining altitude change to achieve the 
commanded y (if no velocity change is called for during the trajectory, this altitude 

change will be the total Ah from the initial to the commanded y) . The values of 

DELHG and DELHV (which are zero if there is no speed change) are summed to compute the 
total altitude change, DHTOT, required to achieve the commanded y. The logic of this 

paragraph is skipped if the y change is completed first; in the latter case DHTOT 

is simply set to DELHG. 

We have now computed che altitude change during the beginning parts of the tra- 
jectory (Ah) from the start of the trajectory to P a (see fig. Bl) . In order to com- 
pute the altitude change during the final portion of the interval, we must first 
compute the velocity at the beginning of the final pitchover maneuver at P 2 in 
figure Bl. If the aircraft completed the velocity change before the y change in the 
initial portion (prior to P 1 ) of the interval, the velocity at P 2 will simply be 
the equivalent airspeed at the end of the interval. Otherwise, V^q for the pitchover 
altitude at P 2 must be computed with an iterative process. The velocity at P 2 is 
initially set to the true airspeed at the end of the flightpath. The procedure for 



determining altitude changes in the previous paragraphs is then executed, resulting 
in an initial estimate for the altitude change from P 2 to the end of the flightpath. 
This altitude change estimate is used to compute a new initial velocity by subtracting 
the altitude-change estimates from the final altitude and determining the velocity at 
the resultant altitude. This velocity is determined by first calculating the altitude 
change needed for the aircraft to complete its required velocity change and then add- 
ing that altitude change to the altitude reached when y e is initially intercepted 
at P x . If the resultant altitude is less than the initial estimated altitude, the 
speed change has been achieved by the time the pitchover maneuver started at point 
P 2 , and no new computation is necessary. If the resultant altitude exceeds the ini- 
tial estimate, the velocity at the initial estimated altitude is computed using the 
formula 


V, = v^2k. Ah/ 1 sin y I + V? 

f 2 1 C 1 X 

and the altitude change computation process is iterated again. The process is con- 
sidered complete when two successive altitude estimates differ by less than 1 ft. 

We now have altitude-change values for both the initial and final portions of 
the flightpath (there is, of course, a constant flightpath angle in the middle portion 
of the flightpath in most cases). If, as with Case 4, the sum of the altitude 
changes for the initial and final parts of the interval total more than the entire 
altitude change desired during the interval, the value for the sin y G is reduced by 
multiplying by the ratio of the total altitude change desired over the sum of the 
initial and final altitude changes, and the resultant y c is used to recompute new 
values for the initial and final altitude changes. 

When the computation of the sum of the initial and final altitude changes is com- 
pleted, the altitude at the beginning of the final pitchover maneuver (ALTGAM) is 
computed by subtracting the final altitude change from the final altitude at the end 
of the interval. Thus, we now have usable values for the commanded value of the 
flightpath angle during the interval and the altitude at which the final pitchover 
maneuver must be initiated. 
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APPENDIX G 


GENERAL DESCRIPTION OF THE QSRA ENERGY CONSERVATIVE GUIDANCE SYSTEM 


The objective of the QSRA ECG system is to synthesize a fly able flightpath from 
the current aircraft states to some desired terminal states. Then, on command from 
the pilot, the system will generate the synthesized flightpath, in real time, and 
provide inputs to an automatic control system or flight director to enable the pilot 
to fly the aircraft along the synthesized path. 

The guidance system is shown in block diagram form in figure Cl, When the pilot 
engages the reference flightpath mode of the autopilot, the guidance executive ini- 
tializes the guidance system. The pilot enters the capture waypoint number to start 
the synthesis computation. The synthesis algorithm first computes the horizontal 
path from the aircraft to the final waypoint and stores the results of this computa- 
tion in the waypoint tables. The horizontal path is affected by the aircraft charac- 
teristics only in the magnitudes of the turning radii that may be used. 

The speed-altitude profile is then found by integrating the equations of motion 
along the horizontal path. This is the most time-consuming part of the synthesis and 
the one in which energy management plays a dominant role. At each point at which a 
change in speed, altitude, or course begins or ends (referred to as a command point), 
a row of data is stored in the command table. An example command table Is given in 
table Cl. Each row in the table contains sufficient information to initialize inte- 
gration in the forward direction to the next command point. Note than IVH = 0 
signifies a null point and the row is skipped. As soon as the synthesis of one com- 
plete flightpath (both speed-altitude and horizontal profiles) is completed, another 
one is started unless the pilot engages the track mode. In the track mode, the 
synthesized reference flightpath is integrated forward in real time from command point 
to command point, and the appropriate information Is generated for the flight director 
and the stability and* control augmentation system (SCAS) . 

All of the operations shown in figure Cl which are inside the dashed lines must 
be done in "foreground" ; that is, these computations have the highest priority and 
must be completed each computation cycle. All other computations, including most of 
the real-time reference flightpath generation, have lower priority and can be a 
background computation, wherein computations are allowed to be interrupted in one 
cycle and completed in subsequent cycles. 

The horizontal flightpath synthesis is illustrated in figure C2. The path con- 
sists of a fixed portion specified by stored waypoints and a capture portion shown by 
the dashed lines. Waypoints are generally used only at the beginnings and ends of 
turns and at the end of the flightpath. Between each pair of waypoints airspeed, 
altitude, or both may change. The speed-altitude profile synthesis causes the speed 
and the altitude changes to terminate at various waypoints. Therefore, it may be 
necessary to insert additional waypoints in order to have the changes terminate where 
desired. 

The horizontal capture flightpath takes the aircraft from its current position 
to one of the fixed waypoints designated by the pilot as the capture waypoint, includ- 
ing touchdown. The complete synthesis of an approach flightpath for the Augmentor 
Wing program requires from 6 to 8 sec (most of it for computing the speed-altitude 
profile), so it was necessary to extrapolate ahead from Px to P 2 (fig. C2) to provide 
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computation time. This extrapolation was done assuming straight, constant-speed 
flight. The pilots, however, continually captured from turns, thus causing large 
initial errors. Therefore, the capture algorithm (subroutines HRZCAP and NEWPSI) was 
modified to allow capture from turns and from non-zero flightpath angles. The algo- 
rithm has also been modified to allow capture paths consisting of three circular arcs 
in addition to the turn-straight-turn paths for improved close-in path generation. 

If the aircraft deviates too far from the reference (as in case of ATC vector) , 
the pilot can disengage and recapture the fixed flightpath at a suitable waypoint. 

(The pilots have expressed an interest in recapturing between waypoints, but this is 
not included in the present implementation.) A number of other options are possible 
for use of the capture algorithm including holding patterns, path stretching, and 
go-around. Their feasibility depends on the computation time for the speed-altitude 
profile and the development of a suitable method for the pilot to control the system. 

The speed-altitude profile is found by integrating the equations of motion for- 
ward from the aircraft along a portion of the flightpath and then backward from the 
final point to the point where forward integration ended. This procedure is explained 
in the text for a simple approach flightpath. Airport-to-airport paths would use 
forward integration to a waypoint near the destination runway and then backward from 
touchdown. This forward /backward procedure ensures the ability to reach the desired 
end conditions, provided the switch-over point is a sufficient distance from the final 
point . 


Synthesis of Horizontal Path 

The fixed portion of the horizontal path is specified by input constants for the 
cartesian coordinates; radius of turn; and maximum, minimum, and nominal equivalent 
airspeeds (ref. 6). The guidance system (subroutine TWOD) starts with the final way- 
point and calculates the course at the preceding waypoint, the angle turned through, 
and the arc length or straight-line distance traversed. Then it iterates in a similar 
manner from waypoint to waypoint until it reaches the capture waypoint. This compu- 
tation is done only once after each selection of a capture waypoint by the pilot. 

The capture path must be continually updated as the aircraft moves until the pilot 
activates the track mode. 

The generation of the capture path (subroutine HRZCAP) is discussed in detail in 
reference 2. The resulting capture path may consist of a circular arc, followed by a 
straight-line segment, followed by a second circular arc, or it may be made up of 
three circular arcs. All of the information about the horizontal path needed for 
subsequent computations Is stored in the waypoint tables for both the capture and 
fixed portions (see table C2) . The first three elements of each table are used for 
the capture flightpath and subsequent ones are used for the fixed path. Only the 
numbers of the fixed waypoints are displayed to the pilot, and these numbers are 
reduced by three from the indices of the waypoint tables. 


Speed-Altitude Profile Synthesis 

An overall description of the technique is given here. The derivations of equa- 
tions, a description of the sources, and manipulation of aero/propulsion data were 
given in appendix A. 
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In order to assure an acceptably flyable flightpath, certain constraints must be 
satisfied. These can be divided into three categories: 

1. Operational constraints: maximum and minimum aerodynamic flightpath angles, 

maximum and minimum airspeed rates, maximum magnitude of normal acceleration, and 
maximum bank angle. These are input constants and may be modified by the pilot 
through the keyboard. 

2. Geometric constraint: maximum magnitude of flightpath angle rate, changing 

from the initial value of y to the commanded y (y c ) during the altitude change 
between two waypoints. If this is not possible because of y c being too large 
(fig. Bl), y c is reduced until a sufficiently small value is found, and the altitude 
at which the final pitchover maneuver begins (hy) is computed and stored. 

3. Aircraft and engine constraints: the result of limited capability of the 

aircraft and engines, including safety constraints. As an example, a flightpath 
angle meeting the constraints in (1) and (2) may not be within the flight envelope of 
the aircraft for a particular flight condition. 

The speed-altitude profile synthesis is illustrated in block diagram form in 
figure C3. The speed-altitude synthesis executive (subroutine VHTSYN) transfers data 
for two adjacent waypoints from the waypoint tables into the synthesis algorithm 
(subroutine INTEG) . In the discussion of the synthesis, these will be referred to as 
the current and next waypoints regardless of the direction of integration. Thus, for 
backward integration, waypoint N is the current waypoint while N-l is the next 
one. 


Closed-Form Equations 

After the initialization in INTEG, a test is made to determine whether there is 
a change in the altitude between the current and next waypoints. If so, starting with 
an initial value of y c , at the maximum operational constraint discussed earlier, 
iteration is used to find a y c for which the change from the initial value of y 
to y and then to the final y can be made during the specified altitude change 
between the waypoints (fig. Bl). The altitude hy at which the pitchover to the 
final y begins is also stored. 

The next step in the synthesis is to integrate the equations of motion, usually 
backward in time, in order to establish the speed, altitude, and heading as functions 
of distance along the reference flightpath. Since the y achieved may be limited in 
subroutine STPINT and not reach the y c computed in the closed-form solutions, the 
pitchover altitude is recomputed shortly before the altitude, computed by integration, 
reaches it. The same closed-form equations are used, but y c is replaced by the 
actual achieved value of y, and the revised value of ALT GAM is stored for use in 
NOMTRJ if the integration is being carried out forward in time. 

If the integration is backward in time, as is the case with most of the synthe- 
sis, the value of ALTGAM obtained from the closed-form equations will not be appro- 
priate for the real-time (forward integration) reference flightpath computed in 
subroutine NOMTRJ. This can be seen from figure Bl. If integration is backward from 
left to right, the pitchover should begin at P 2: , while for forward integration (from 
right to left), the pitchover altitude is that at P x . Again, the aircraft may not 
actually be able to achieve y c , so when the backward integration reaches the alti- 
tude computed for ? 1 by the closed-form solution, the closed-form computation is 
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repeated using the value of y actually achieved. This value is stored for use in 
the real-time computation in subroutine NOMTRJ. The resulting y c is only a geomet- 
rical limit on |yj, however, and the aircraft may be incapable of achieving it 
because of aircraft and/or engine constraints. For this reason, as hy is approached 
during the integration, it is recomputed using the current value of y. There is a 
separate part of subroutine INTEG which computes y c and hy for this latter case. 

The derivatives of the closed-form equations and a further discussion of their 
use are given in appendix B. 


Real-Time Reference Generation 

The real-time reference generation is illustrated in block diagram form in 
figure C4. The algorithm is divided into two phases. The first (subroutine NOMTRJ) 
essentially repeats the fast-time integration (but not the closed-form calculations) 
done in the synthesis. In this case, however, the integration is all done forward in 
time. The integration between successive command points is initialized from the 
command table. Large integration step sizes (1 or 2 sec) are used as in the synthe- 
sis. In the flight program the integration is done in background. The resulting 
data are stored in the 15 x 6 array, RDERIV. Each column of RDERIV corresponds to 
the end of one of the long integration steps. When a command point is reached, the 
step size is adjusted to the distance to the command point. At the command point the 
speed, altitude, and heading are set to the values stored in the command table. 

The final step in the real-time reference generation is the integration in real 
time between the columns of RDERIV (entry N0MTR2) . Distance is used as the indepen- 
dent variable, and in the flight program, the step size is set to the projection of 
the distance traveled by the aircraft in one frame time (currently 0.1 sec) projected 
on the reference path. When the distance corresponding to one of the columns of 
RDERIV is passed, the step size is reduced to the exact value needed and the aircraft 
states are corrected to the values stored in RDERIV. The remainder of the small 
integration step is added to the one immediately following. The horizontal cartesian 
coordinates are generated only in this final real-time computation, and at waypoints 
they are corrected to the values in the waypoint tables generated by the horizontal- 
path synthesis. 

At some command points the synthesis may call for a step change in V a /g, which 
may in turn require step changes in flaps and, if | sin y + V a /g| exceeds the limits 
on E n , in y a . These step changes are smoothed out by computing perturbations in 
speed, altitude, flaps, and y a to be added to the values from the tables for a lead 
distance on each side of the command point (see fig. C5) . The perturbations are com- 
puted by assuming linear changes in y a , V a /g, and flaps over the distances needed, 
such that for each variable, half of the step change is achieved at the command point. 


The Carrier Landing Modifications 

The program has been modified to account for straight, constant-speed motion of 
the aircraft carrier. This is accomplished by adding a velocity which is equal in 
magnitude and opposite in direction to the ship velocity to the estimated winds. This 
is done in subroutine WINDAT for all the tabulated values of estimated winds. 

This procedure results in the aircraft motion in the basic program being defined 
relative to the moving, ship-fixed coordinate system. As a result, the inertial 
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flightpath angle y*j- is relative to the ship-fixed system. However, it appears 
desirable to limit the flightpath angle, YpE 5 w i t ^ L respect to the Earth-fixed inertial 
system. If we define Vjj? as the velocity of the aircraft with respect to the 
Earth-fixed system, Vj as the velocity with respect to the ship-fixed system, and 
V a as the true airspeed vector , then 

V IE sin y ie = V ][ sin Yj = sin Y a (Cl) 


where y a is the aerodynamic flightpath angle. From equation (Cl), 


XE 

sin Y a = y ~ sln Y ie (C2) 

v a 

Equation (C2) is used to compute the limits on sin y a from the limits on sin YqE» 
which are input constants. 

Now define the velocity vectors. The first set are new definitions of variables 
in the basic program. 

V true airspeed vector (no change) 
a 

V velocity of the wind with respect to the ship (used to be with respect to Earth) 

Vj velocity of aircraft with respect to ship (used to be with respect to Earth) 

New variables are: 

V velocity of ship with respect to Earth 
SH 

V „ velocity of wind with respect to Earth (input estimated values) 
wE 

V velocity of aircraft with respect to Earth 
The following are relationships beweeen the vectors 


or 



(C3) 


v = V + v 
IE I SH 


(C4) 


From eqation (C4) 

V IE " V V I + 2V I V SH C0S(HR) + V SH <C5) 

where HR is the course of the aircraft with respect to the ship. 
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TABLE 1.- WAYPOINTS FOR PATH NO. 1 


Waypoint 

No. 

EPSLN 

X, 

ft 

Y, 

ft 

D 

Turn 

radius, 

ft 

Nominal 

airspeed, 

kt 

Maximum 

airspeed, 

kt 

i 

1 

-59311 

52166 

-2000 

0 

140 

160 

2 

1 

-57030 


-2000 

0 

140 

160 

3 

1 

-52864 

44869 

-2000 

9258 

140 

160 

4 

1 

-41767 

12440 

-1200 

0 

140 

160 

5 

1 

-35745 

7164 

-1200 

8033 

140 

160 

6 

I 

-26127 

5236 

-1200 

0 

80 

90 

7 

1 

-17872 

3582 

-1200 

0 

65 

75 

8 

1 

0 

0 

-60 

0 

65 

75 


TABLE 2.- WAYPOINTS FOR PATH NO. 2 


Waypoint 

No. 

EPSLN 

X, 

ft 

H 

n 

Turn 
radius , 
ft 

Nominal 

airspeed, 

kt 

Maximum 

airspeed, 

kt 

i 

1 

-59311 

52166 

-2000 

0 

140 

— 

2 

1 

-57030 

50160 

-2000 

0 

140 


3 

1 

-52864 

44869 

-2000 

9258 

140 

160 

4 

1 

-41767 

12440 

-1200 

0 

140 

160 

5 

1 

-35745 

7164 

-1200 

8033 

140 

160 

6 

1 

-17872 

3582 

-1200 

0 

65 

75 

7 

1 

0 

0 

-60 

0 

65 

75 


TABLE 3.- WAYPOINTS FOR PATH NO. 3 


Waypoint 

No. 

EPSLN 

X, 

ft 

Y, 

ft 

z, 

ft 

Turn 
radius , 
ft 

Nominal 

airspeed, 

kt 

Maximum 

airspeed, 

kt 

i 

1 

-17872 

3582 


0 

65 

75 

2 

1 

0 

0 

■S3 

0 

65 

75 



TABLE 4 

.- WAYPOINTS FOR PATH NO. 4 


' '■ ■ ' 

Waypoint 

No. 

EPSLN 

X, 

ft 

Y, 

ft 

z, 

ft 

Turn 
radius , 
ft 

Nominal 

airspeed, 

kt 

Maximum 

airspeed, 

kt 

i 

i 

-5842 

1171 

-919 

0 

65 

75 

2 

i 

0 

0 

-60 

0 

65 

75 
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TABLE 5.- WAYPOINTS FOR PATH NO. 5 


Waypoint 

No. 

EPSLN 

X, 

ft 

Y, 

ft 

z, 

ft 

Turn 
radius , 
ft 

Nominal 

airspeed, 

kt 

Maximum 

airspeed, 

kt 

1 

2 

1 

1 

-2974 

0 

596 

0 

-919 

-60 

0 

0 

65 

65 

75 

75 


TABLE 6.- FLIGHTPATH RESULTS SUMMARY 


■ 

1 

Distance 
flown, ft 

Flight 
time, sec 

Fuel 

used , lb 

1 

95,062 

747 

1148 

2 

95,062 

701 

1080 

3 

88,254 

676 

1070 

4 

85,999 

529 

685 

5 

96,314 

442 

552 


34 


i 






























TABLE Al.- ENGINE PARAMETERS 


Index 

Gross 

thrust 

Ambient 

pressure 

ratio 

constraint 

Aircraft 

weight 

i 

Maximum 

1.000 

Maximum 

Minimum 

2 

Maximum 

.688 

Maximum 

Minimum 

3 

Maximum 

.459 

Maximum 

Minimum 

4 

Minimum 

1.0 

Maximum 

Maximum 

5 

Minimum 

.688 

Maximum 

Maximum 

6 

Minimum 

1.0 

Minimum 

Maximum 

7 

Minimum 

.688 

Minimum 

Maximum 


TABLE A2 LIMITS ON GROSS THRUST DUE TO 
LIMITS ON ENGINE CORE SPEED 


Mach no. 

t gmax^ 5 amb’ lb 

t gmin /5 amb’ lb 

0 

25,400 

8,200 

.1 

25,400 

8,200 

.2 

27,000 

9 , 100 

.3 

29,520 

10,900 


TABLE A3.- VERSUS N //¥~ AND 

G AMB 12 

MACH NO. (ONE ENGINE) 


Nj/Ze^ 

o 

II 

S 

M = 0,1 

M = 0.2 

M = 0.3 

80 

4416 

4539 

4824 

5235 

82 

4700 

4823 

5094 

5507 

84 

4984 

5108 

5364 

5839 

86 

5268 

5392 

5659 

6172 

88 

5559 

5697 

5981 

6504 

90 

5878 

6014 

6303 

6837 

92 

6197 

6332 

6625 

7167 

94 

6516 

6650 

6947 

7499 









TABLE Cl.- SAMPLE COMMAND TABLE WITH COMPUTER- VARIABLE DEFINITIONS 


IWPT 

IVH 

KCOR 

T 

D 

GAMMA 

VADTG 

FLAPS 

PHID 

ALTS 

VIASS 

DLPHI 

DLFAC 

1 

5 

1 

298.9 

50000.0 

0. 

00 

0. 

00 

0 

.00 

0.00 

500 

.00 

236 

.60 

47794 

48252 

1 

0 

4 

288.9 

47793.9 

0. 

00 



0 

.00 


0.00 

0 

.00 

47794 

48252 

1 

0 

0 

288.9 

47793.9 

0. 

00 



0 

.00 


0.00 

0 

.00 

47794 

48252 

1 

0 

0 

288.9 

47793.9 

-0. 

04 



0 

.00 


0 

.00 

0 

.00 

47794 

48252 

2 

2 

0 

288.9 

47793.9 

0.03 



-59.00 


500 

.00 

236 

.60 

39696 

40797 

2 

0 

1 

249.1 

39695.6 

0. 

00 



0 

.00 


500 

.00 

202 

.80 

39696 

39696 

2 

5 

1 

249.1 

39695.6 

0. 

00 



0 

.00 


500.00 

202 

.80 

14739 

14739 

2 

1 

1 

115.4 

14738.6 

0.01 

-0. 

04 

65.97 


500 

.00 

202 

.80 

9516 

10552 

3 

0 

0 

84.6 

9516.0 

0. 

00 

0 

00 

0 

.00 


0 

.00 

0.00 



9516 

3 

0 

0 




















3 

0 

0 






. 














3 

0 

0 




















8 

0 

0 




















8 

0 

0 





1 

f 




f 


1 

f 

\ 

r 


f 


f 

8 

3 

0 


f 

1 

f 

-0.05 


t 

0.66 


996 

.76 

169.00 

8405 

8705 

8 

1 

1 

77.5 

8404.8 

0 

.05 

-0. 

.03 

-12 

.41 


917.17 

169.00 

3800 

3919 

9 

0 

0 

40.6 

3800.0 

0 

.00 

0, 

.00 

0 .00 


o 

.00 

0.00 

3800 

3800 

9 

0 

0 

40.6 

3800.0 

0 

.00 

0, 

.00 

0.00 


0.00 

0.00 

3800 

3800 

9 

0 

0 

40.6 

3800.0 

0.00 

0, 

.00 

0.00 


0 

.00 

0 

.00 

3800 

3800 

9 

3 

0 

40.6 

3800.0 

0 

.00 

0, 

.00 

0.00 


600 

.35 

111.54 


0 


0 


Variable definitions: 

IWPT Waypoint number 

IVH Defined quantity (in INTEG) showing type of flight path change: 

0 - null point, no change in state 

1 - speed and altitude changing 

2 - speed only changing 

3 - altitude only changing 

4 - constant speed level turn 

5 - constant speed, straight and level flight 

KCOR Flag used to correct integrated path to correspond to horizontal capture 
path 

T Time to go to touchdown 

D Distance to touchdown 

GAMMA Aerodynamic flight path angle 

VADTG Airspeed rate in g T s 

FLAPS Sum of outboard and USB flap settings 

PHID Bank angle at command point 
ALTS Altitude at command point 

VIASS Indicated airspeed at command point 

DLPHI Distance at which roll command changes (includes lead distance) 

DLFAC Distance at which flaps and V a /g change (includes lead distance) 


j 
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TABLE C2.~ SAMPLE WAYPOINT TABLE WITH COMPUTER-VARIABLE DEFINITIONS 
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A START OF CAPTURE PATH 
B END OF FIRST CAPTURE TURN 
C START OF SECOND CAPTURE TURN 
D CAPTURE WAYPOINT 

Figure 1.- Horizontal profile of capture trajectory. 
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X* 






X, ft 



Figure 6.- Path No 


5 horizontal profile with waypoints. 
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ROLL PITCH FLIGHT PATH 

ATTITUDE, HEADING, ATTITUDE, ANGLE, ALTITUDE, 




DISTANCE TO TOUCHDOWN, ft 

Figure 7.- Synthesized data for Path No. 1. 
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OUTBOARD ENGINE INDICATED 
FUEL USED, FUEL RATE, USB FLAPS, FLAPS, FAN SPEED, AIRSPEED, 
lb Ib/sec deg deg percent knots 
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Figure 7.- Concluded, 
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OUTBOARD ENGINE INDICATED 
FUEL USED, FUEL RATE, USB FLAPS, FLAPS, FAN SPEED, AIRSPEED, 
lb Ib/sec deg deg percent knots 



20,000 40,000 60,000 80,000 100,000 

DISTANCE TO TOUCHDOWN, ft 

Figure 8.- Concluded. 
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OUTBOARD ENGINE INDICATED 
FUEL USED, FUEL RATE, USB FLAPS, FLAPS, FAN SPEED, AIRSPEED, 
lb Ib/sec deg deg percent knots 




800 
400 

0 20,000 40,000 60,000 80,000 100,000 

DISTANCE TO TOUCHDOWN, ft 



Figure 10.- Concluded. 
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INDICATED INDICATED 

FUEL RATE, ALTITUDE, AIRSPEED, FUEL RATE, ALTITUDE, AIRSPEED, 
lb/sec ft knots Ib/sec ft knots 
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0 1 2 3 4 5 6 7 

C L 

(a) Total energy rate (E^) versus lift coefficient (C L ) . 
Figure A3.- Constant thrust coefficient (C^,) . 
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(c) Flap setting (fi™--.-,) versus lift coefficient (C T ) . 

rK.hr Xj 


Figure A3.- Concluded. 
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(a) Total energy rate (E T ) versus lift coefficient (C^) 
Figure A4.- Constant thrust coefficient (C T ) . 







• 

C T 

MAX 

5 AMB 

1.000 

ALT, ft 
SEA LEVEL 

• 

E-r 

o 

MAX 

0.687 

10,000 

, 1 MAX 

Et 

□ 

MAX 

0.459 

20,000 

t min 

A 

MIN 

1.000 

SEA LEVEL 


V 

MIN 

0.687 

10,000 




Figure A4.- Continued. 




C L 


(d) Thrust coefficient (C T ) versus lift coefficient (C^) . 
Figure A4.- Concluded. 
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Figure A6.- Corrected engine fan speed (N x / /0^") for exhaust gas temperature (MGT) 

of 920° C versus temperature change (AT). 


7 C TOO LARGE 


PATH IF A/C CAPTURES 
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Figure Bl.- Closed-form solutions illustration. 




i 


67 



CALLED BY AUTOPILOT EXECUTIVE 



L 


J 


Figure Cl.- QSRA ECG guidance system. 
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Figure C2.- Horizontal flightpaths displayed on cockpit map display. 
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Figure C3.- Speed-altitude profile synthesis. 
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Figure (J4.- Real-time reference generation. 
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Figure C5.- Perturbations in V a /g, 7 5> and flaps resulting from step changes 

in V a /g, 
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